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We study the dewetting of a porous plate withdrawn from a bath of fluid. The micro- 
scopic contact angle is fixed to zero and the flow is assumed to be parallel to the plate 
(lubrication approximation) . The ordinary differential equation involving the position of 
the water surface is analysed in phase space by means of numerical integration. We show 
the existence of a critical value of the capillary number rjU /"/, above which no stationary 
contact line can exist. An analytical model, based on asymptotic matching is developed, 
that reproduces the dependence of the critical capillary number on the angle of the plate 
with respect to the horizontal for large control parameters (3/2 power law). 



1. Introduction 

When sea retreats from the shore, sand structures appear as solid granular particles are 
transported via the liquid. Liquid motion and particularly film retraction on an erodible 
medium are known to create impressive erosion patterns, s uch as sand ripple s for oscilla- 
tory waves JStegner fc Wesfrei dl ll999HScherer et. raZjITflflfo or sand furrows l|Paerr et all 



2003; Schorghofer et al. 2004). The case of liquid retraction from a granular bed can be 
understood as a dewetting dynamics on a porous erodible bed. Such physical phenomena 
have been repro duced in the labor atory by pulling a plate covered with a bed of grains out 
of a liquid tank ()Paerr et al.l2003|) . This situation is similar t o the well-known experim ent 
investigating a moving contact line on a non-porous plate l|Blake fc Ruschaklll979(l . In 
this latter case, a contact line exists for small removal speed U, whereas for higher speed 
(above a well-defined critical value U cr ) a macroscop ic water film (the so-calle d Landau- 
Levich -Der i aguin film (denote d LLP later on ), see lLandau fc Levichl l|l942j) : IPeriagfnl 
( 1943)) covers the whole plane l)Eggersll2004al) . We propose here to investigate this tran- 
sition for a saturated porous medium, in connection with recent experiments involving 
granular materials bv lPaerr et all l(2003h . There, a motor-driven plane, covered with a 
granular layer, is withdrawn from a water tank at constant speed U. The solid plane is 
tilted to an angle 9. At high enough velocity, erosion river networks and mudflows are 
observed, whereas only light patterns appear at smaller speed. We investigate the loss 
of a static contact line and seek to relate it to the transition between various erosion 
regimes. We therefore seek the critical velocity above which no static contact line can 
exist on a granular bed. Below this critical velocity, almost no grain motion is observed 
so that we identify the granular bed with a rigid porous medium. Pewetting on a porous 
mediu m has already b e en studied in different configurat i ons (s ee [Raphael fc de Genned 
l|l999j) : lAradian etall i|200(]j) ; iBacri fc Brochard-Wvartl l|200l[) and references herein) 
The c ase of a porous plate removed from a liquid was studied bv lRaohael fc de Gennesl 
l)l999l) . but the focus was on the spatial evolution of the LLP film, the existence of which 
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Figure 1 . A porous plate of conductivity k is being withdrawn from a liquid bath with speed 

U at angle 6. 



was assumed. The contact line dynamics was also studied IjAradian et all 2000) for a wet 
horizontal support in which the liquid is sucked in the dry porous medium. 

From a more theoretical point of view, the problem of a moving contact line on a 
porous solid is pointed out bv lde Penned lll985h as a natural regularization for the contact 
line dynamics equations. The bulk liquid flow through the porous solid indeed removes 
the usual stress singularity that one would encounter at a contact line with a no-slip 
condition l|Dussan fc DavisHl974|) . No additional assumption, such as the introduction 
of a Navier slip at microscopic scale, is then required. However, another question arises 
when considering a porous medium: what is the relevant condition for the contact angle 
at the contact line? As discussed below, we propose here for a saturated system to take 
a zero contact angle. 

The paper is organized as follows: in the next Section, we use the lubrication approx- 
imation to deduce the equation for the static interface shape, both for the contact line 
and a zero-flux LLD film. In Section |21 we exhibit the transition between these two con- 
figurations as the pulling velocity increases by means of a shooting method. Then, we 
propose to interpret the solutions in the framework of dynamical systems (Section 0J. 

2. Principles 

2.1. Lubrication approximation 

Our approach seeks to determine the velocity (if any) above which the static contact 
line can no longer exist in a granular bed withdrawal experiment. Below this velocity, 
we can consider that the grains almost do not move relatively to the withdrawn plate. 
Thus the granular material is represented by a non-erodible porous medium (see Figure 
^1 of permeability fc, and we only have to investigate the stationary problem. The fluid 
is characterized by its density p, dynamical viscosity rj and surface tension 7. Assuming 
invariance in the ^-direction, we consider the two-dimensional problem where the water 
surface is described by the function h(x). For the plate velocities pertinent to the problem 
(typically 0.5 cm. s -1 ) and the estimated porosity of the granular bed (k ss 10~ 12 m 2 ) we 
can consider that the porous medium remains fully saturated with water at any distance 
from the free water level. 

We will restrict our analysis to small angles 9 so that the lubrication approximation 
can be employed {i.e. 9 , \\h'\\ <C 1 where the ' stands for the x-derivative). Only the 
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IE-component u of the velocity has to be taken into account, for which Poiscuillc profile 
is assumed, with a vanishing tangential stress on the gas side: 



du 
dy 



0. 



(2.1) 



y=h(x) 



Another boundary condition has to be written at the porous surface. The classical no- 
slip condition, as required at the solid-fluid interface on an impermeable plate, leads to 
the following equation for h(x) (its derivation is similar to the one presented in Appendix 
IS): 



ti 



3Ca 



(2.2) 



where lengths h and x have been made dimensionless by the capillary length l c 
y/'y/pg (l c ~ 2.8 mm for water) and Ca is the capillary number, defined by 



Ca: 



7 ' 



Equ ation 112.211 and any der ivative of its solutions are singular at the cont act line, where 
h = (Du ffy fc Wilsonlll997|) . For a non-porous surface, a short-l ength regularization i s 
invoked coming either from effective slip near the contact line ( Huh fc Scrivenll l97li). 
the existence of a pre-wet ting liquid film and Van der Waals forces l | de Gennesl 1984; 
iHervet fc de Genneslll984Tl or a " diffuse interface model" l|Seppecherlll99fil) . Such a regu- 
larization always involves a microscopic cut-off length (on the order of 1 nm) below which 
it is claimed that hydrodynamics fails. The Navier slip condition i s then mostly used i n 
numerical simulations investigating moving contact line problems ijRenardv et a/.ll200l|) . 
This condition reads at the solid-fluid interface: u — U = A^du/dy at y = where Ajv 
is the cut-off length. If Ajv is the rescaled cut-off length (that is A^v = Ajv/Z c ), Equation 
12.21) becomes 



b!" -h' + 9 = 



Ca 



h 2 /3 + X N h' 



(2.3) 



Equation 12.. 'ill can be numerically solved and analytically approached. A contact line 
is then found to exist as long as the capillary number is smalle r than a critic al value, 
above which a macroscopic LLD film is deposited on the solid (|Eggersll2004l) . Notice 
however that not all the singularities discussed above are suppressed by the Navier slip 
condition since the capillary pressure still diverges at the contact line (see Appendix 0. 



2.2. The case of porous solid 

A porous solid allows for both interfacial slip (first proposed bv lBeavers fc Joseph! l|l967f) ) 
and bulk flow. Using the Brin kman equation to describe the flow inside the porous 
medium. iNeale fc Naderl l|l974|) showed that, for a homogeneous porous media, the mag- 
nitude of the slip is proportional to the prevailing shear stress: 



l»=0 



L P\y=0 



y/k du 
a dy 



(2.4) 



y=o 
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where k is the permeability of the solid, a a coefficient of order one, and u p the velocity 
of the fluid in the porous medium. Darcy's law holds in the solid so that 

u P + U = --(^-pg6), (2.5) 



r\ \dx 

and finally the Equation hereafter describes the shape of the steady fluid film under 
withdrawal (the detailed derivation is presented in Appendix lB|l : 

ti" -h' + 9= u9/o ° a . -, (2.6) 

h 2 /S + Xh/a + A 2 v ; 

where A = \fkjl c . Such an equation is similar to those studied for the contact line on 
a solid surface , using specific boundary conditions at the solid surface i|Eggersl 12004 al 
iHockindEloTI) . It has been shown (in the case of a plate pushed into water) that the 
details of the regulariz ation do not influence the far- field fluid flow as long as the cut-off 
length is small enough l|Eggersll20046|) . However, an importan t difference in our c ase lies 
in the typical values of A involved in porous media (s» 10~ 2 in lDaerr et al\ l)2003|) ) to be 
compared with 10~ 6 for regular solids. 

Here we would like to point out that in recent papers ijMaurer et a?j20o3lHadiiconstantinoul 
120031) a second order-slip law is used to model the flow of gases at large Knudsen numbers. 
This boundary condition (adapted to the present notations) reads 



Ay^ + U = k C — 



. d 2 u 

- Ctc A C , 
y=o °V y=0 

where etc is a positive coefficient of order one, and Ac? is a slip length of the same 
order than the mean free path of the gas. If such a boundary condition were used in the 
case of a contact line, again Equation l|2.(jl) would be obtained. 

2.3. Boundary conditions 

The limit for large positive values of x is well-defined: the water surface is horizontal far 
from the plane, that is 

h{x) ~ Ox. (2.7) 

x — >oo 

For the two remaining boundary conditions, two different cases will be studied, de- 
pending on whether a contact line is formed between the water surface and the solid 
plate, or if a film of water remains on the solid surface. In the first case, the water level 
vanishes at the origin and a contact angle #o is usually imposed; the contact line set of 
boundary conditions is 

h'(o) = e . (2 - 8) 

Following iRaphael fc de Gennesl l|l999|) . we will hereafter consider that the contact 
angle O is zero for dewetting on a porous media. We argue indeed that for saturated 
porous media, the liquid film wets completely the surface, leading to an effective zero 
contact angle. 

When LLD film starts at the meniscus, the only boundary condition known a priori is 



lim h(x) = hf 



(2.9) 
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where ft,/ is a constant solution of Equation H2.6J) . Notice that Equation (|2.6|l stands 
only for zero- flux films, that neither add nor withdraw water from the tank. We omit 
the important case wher e a LL P film is continuously growing with time, a case that was 
investigated bv lHockind (l20Q]J) . Our study is therefore relevant to determine the loss of 
a static contact line solution, without any information about the dynamics. The stability 
of the solution as well as the time-dependent dynamics of a moving meniscus cannot be 
studied at this stage and will be the purpose of further work. However, we will see in 
Section 0| that the film solutions of (|2.6|) satisfying (|2.9I) play an important role in the 
dynamical system describing our solutions. Finally, for a solution to be acceptable, the 
water level must always lay above the porous medium: Vx, h(x) > 0. 

2.4. Parameters 

The parameter a comes fr om the detailed modelisation of the interf ace slip flow jBeaver s fc Joseph! 
Il967t iGoveau et alJl2003IL and varies generally between 0.1 and 4 IjNeale fc Naderl ll9741. 
Notice for instance that if a is smaller than Equation 1|2.6[) becomes singular for 

some negative values of h. However, for the sake of simplicity the coefficient a is fixed 
to one in the present study. Now, for a = 1, Equation (|2.6I) is an ordinary differential 
equation (ODE) with three parameters: 8, Ca and A. In fact, one should notice that this 
equation is only a two-parameter ODE: defining h* = h/X, Equation (|2.(il) becomes 

h '"- k > +l> ' = ky3%, + v < 2 ' 10 > 

where 9* and Ca* are defined as follow: 

e* = e/x = 



Ca* = Ca/A 3 = 



7 
pgk' 



( w fc) 3/2 ' 

Equation (|2.1U|) is the one we will study later on, but we will omit the * on h* for 
the sake of readability. These nondimensional parameters were choosen because they are 
proportional to the two experimental parameters that may be easily and continuously 
tuned, namely U and 9. If the permeability of the porous solid is extremely low, both 
9* and Ca* tend to infinity, as well as the ratio C&*/8* = Ur)/(6pgk). In this case, the 
velocity inside the po r ous m edium is extremely slow compared to U. For the experimental 
study of lDaerr et all (|2003|) . the rescaled capillary number Ca* is of the order of 10 5 . 

2.5. Hydrostatic solutions 
Any solution which respects the boundary condition l|2.7|l for large x verifies 

lim h(x) = +oo, 

x — >+oo 

thus for large x, Equation (|2.1U|) becomes 

ti" - h! + 6* = 0. (2.11) 

The behaviour of the water surface at large x, hereafter denoted by ftoo, is directly 
obtained from Equation l|2.11|l : 

h x {x) =A OD + 8*x + (9* - 0* p ) exp (-x), (2.12) 
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Figure 2. Numerical solutions of Equation l|2.1()^ for different capillary numbers. The rescaled 

tilt angle 8* is fixed to 1 



where and 9* p are two constants, corresponding respectively to the length of the dy- 
namical meniscus and to the so-called apparent contact angle (note that the dimensional 
apparent contact angle is actually 0* p X). Even though the lubrication approximation is 
not expected to hold for large x since the water level is not small anymore, Equation 
H2.11f> leads to the classical static meniscus solution (remember that x has been scaled 
by the capillary length l c ). Consequently, we may consider that Equation p.lOfl holds at 
any position on the s-axis. 



3. Numerical results 

3.1. Contact line solutions 

To seek steady contact line solutions, Equation l|2.10|l may be solved numerically using a 
finite-difference algorithm. In the case of a contact line, two boundary conditions may be 
fixed at x — by the contact line conditions 1)2.8(1 . The third condition comes from the 
flat water level at infinity l|2. 7JI . With the new notations, we end up with the following 
system: 

ti" -h' = f(h), 

h(0) = h'{0)=0, (3.1) 
h(x) ~ 9*x, 

X — >oc 

where 

We use a shooting method (see lMannevillell)l99nln . varying the initial curvature h"(0) 
in order to find the numerical solution which corresponds to the hydrostatic condition at 
large x. Some numerical contact line solutions to Equation (|3.1|l are presented on Figure 
121 for different Ca* at fixed 6*. 
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Figure 3. Curvature at the origin versus rescaled capillary number, for a contact line solution of 
Equation l2.10l The rescaled tilt angle 9* is fixed to 1. No solution is found for capillary numbers 
higher than Ca* w 1.0247. 



As shown on Figure |2 the contact line zone is somehow streched as the capillary num- 
ber is increased. In other words, the curvature h"(0) at the origin tends to zero as Ca* 
tends to a critical value Ca*. Above this critical value, no contact line solution can be 
found by the shooting mehod. This transition will be clarified below using the dynam- 
ical system associated to 113.10 . The disappearance of the contact line solution may be 
represented in a kind of bifurcation diagram, plotting the curvature at the origin against 
the capillary number, as shown in Figure EJ Notice that, even though the curvature at 
the origin tends to zero as Ca* is approached, the contact line solution does not become 
unrealistic for Ca > Ca* owing to a negative initial curvature, but rather disappears 
by a bifurcation. Above the critical capillary number, no matching exists between the 
behavior of the solution at the contact line and the gravity-capillary solution. 

3.2. Film solutions 

A LLD solution can exist whenever there is a positive value hf such that f(hf) = 0, 
which occurs as soon as Ca* ^ 6*. Two major limitations have to be pointed out for these 
film solutions: first, we might not be able to match this film solution to the hydrostatic 
region with h remaining positive everywhere. Moreover we restrict our analysis here to an 
already-established fi lm of zero mas s flux, whereas transitory and/or finite flux solutions 
should be considered l|Hockina200l|) . and are likely to exist for smaller capillary numbers. 
Numerically, above the critical capillary number Ca*, we have always been able to find 
a zero flux LLD film of thickness hf in the limit x — > — oo that could match to the 
hydrostatic solution without crossing h = (see Figure EJ- Such a solution may be 
numerically approached, using a special shooting method described in Section B.l.ll We 
observed that as the capillary number is decreased, the film surface is shifted down along 
the y-axis (see Figure |HJ), and we may define a second critical capillary number Ca* 2 , 
bellow which the film solution becomes negative in some region. Consequently, if Ca* 2 
is smaller than Ca*, hysteresis may occurs, that is, two solutions, a contact line one and 
LLD film one, may co-exist for same tilt angle and capillary number. 
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Figure 4. Numerical solutions of Equation I2.10H for various capillary numbers, represented 
in the phase space (projected on the (h,h")~ and (h! , /i")-planes). The rescaled tilt angle 6* is 
fixed to 1. The trajectories correspond to the physical solutions shown in Figure|5| for subcritical 
capillary numbers (the three dashed curves), the trajectory starts at a point on the /i"-axis which 
corresponds to the contact line. On the opposite, the solid line correspond to a film solution, 
and thus does not cross the /i"-axis. The insert shows the projection of the film solution in the 
(h! , /i")-plane, at smaller scale. 



4. Dynamical systems interpretation 

4.1. Phase space 

In what follow, we inte rpret and develop the preceding results using dynamical system 
theory l)Strogat7l ll994^. Let us consider the phase space V corresponding to Equation 
l|2.1(J[l . that is, M. 3 with coordinates (h,h',h"). Any solution of 12.10fl is a trajectory of 
V, parametrized by x, which satisfies 

/ h' \ 

X' = F(X) = h" . (4.1) 

V h'+f(h) J 

Some trajectories (the same as on Figure |2J are represented on Figure 0] Notice that 
the film solution (solid curve) winds exponentially around a fixed point on the h-axis. 

4.1.1. Hydrostatic solutions in the phase space 

For large x, following the reasoning of Section l2~5l Equation (|4.1|l becomes linear: 



X'= 1 X- . (4.2) 




Any solution of (|4.2() which satisfies the boundary condition l|2.7[l is included in a 
plane called i?oo, which may be parameterized by x and the apparent contact angle 9* p 
introduced in Section 12.51 Eoo is defined by the equation h' — h" = 9* . 

The solutions of the full Equation l|4.1|) which satisfy l|2.7|l are included in a two- 
dimensional manifold, called W. This manifold tends to for large h. This allows us 
to approximate numerically the LLD film trajectories, for which we impose boundary 
conditions at x — > — oo and x — > +oo. We may indeed use a shooting method with initial 
conditions varying along a constant (large) h line on E^. The boundary condition at 
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x — > +00 is then approximatively satisfied at any step. The shooting method provides 
an approximation of the only solution that remains constant as x tends to —00. 

Figure represents the intersection of W with the (h, h") plane defined by h' = 0, 
obtained by the shooting method, for two different capillary numbers, above and below 
Ca*, for 9* = 2. We observe numerically that the major effect of an increase in Ca* is a 
translation in the higher h direction. The disappearance of the contact line solution may 
be described in the following way: any contact line trajectory is embedded in W, and the 
boundary conditions impose that it starts on the /i"-axis, consequently, it can exist only 
if there is an intersection between W and the /i"-axis. Since the main effect of an increase 
of Ca* on W is a translation along the ft,-axis, this intersection disappears above some 
value Ca* of the capillary number. Thus, the existence of a fixed point creates a separatrix 
on the boundary W, that would otherwise be defined over the whole (h, /i')-plane, which 
allows for the sudden disappearance of its intersection with the /i"-axis. 

4.1.2. Fixed points 

A fixed point X/ in phase space corresponds physicaly to a film of constant height hf. 




The existence and values of fixed points depend on the parameters 9* and Ca*, as 
presented in Tabled I n the following, we will focus on the largest fixed point , since 
it is the only one that may be acceptable physically (that is ft./ > 0). Let us linearize 
Equation (|4.1() around X^: 

X' = J+ (x-x+), 

where Jt is the jacobian of F evaluated at Xt, that is 
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I Condition I Fixed points h 



Ca* < 



0* 



Ca* 



Ca* > 



-1 ± 



1 / 4Ca 
V3 



Table 1 . Existence and values of the fixed points of equation l|4.1|l . 



| Condition ) 


Eigenvalues | 


| Ca* 6 T(6»*) | 


a, b, c | 


|Ca* g <9T((9*)| 


a, -a/2, -a/2 | 


| Ca* £T(0*) \a, 


-a/2 + iQ,, -a/2- iQ \ 



Table 2. Eigenvalues of the jacobian Jj at the largest fixed point, a, b, c and Q are real numbers 

and T(6>*) = (l - y^l - gjW) ; ^± (l + y^l - 55W)] • When the eigenvalues are real, 

they satisfy: a < < b < c. 




J+= 
V f'(hp 

The local behaviour of solutions around the fixed point depends on the eigenvalues of 
J j , which are presented in Table [21 If the eigenvalues are real numbers, one is negative 
and the two others positive. So there is an unstable manifold of dimension two where 
the trajectories tend monotonically to the fixed point as x tends to -co and a stable 
manifold S (separatrix) of dimension one. On the other hand, when the two eigenvalues 
are complex conjugate their common real part is always positive, and the trajectories in 
the corresponding unstable manifold wind around the fixed point while diverging from it 
at exponential rate. In physical space, the fluid surface forms damped stationary waves 
along the plate (see Figure EJ. Therefore, the behaviour of W in the vicinity of the fixed 
point may also be described by this linear expansion. Depending on the parameters, W 
may either be defined over the whole (/i, /i')-plane, or tend to the separatrix S (which 
ends on the fixed point). In the latter case, W winds around the separatrix (as shown 
in Figure EJ or tends to it monotonically. These various regimes are represented in the 
diagram of Figure 
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4.1.3. Critical capillary number 

At any point on the /i-axis, X' is parallel to the ft,"-axis (see equation 1)4. ljl ). Con- 
sequently, the intersection of W with the (h, h") plane (represented on Figure [SJ has 
vertical tangent vectors whenever it crosses the h-axis. This explains the behaviour of 
h"(0) close to the critical capillary number (see Figure |3J), which may be interpreted as a 
saddle-node bifurcation. This property is useful for the numerical determination of Ca* 
at a given 8* : since we know that the second derivative h" (0) must vanish at the critical 
capillary number, we may approximate Ca* by a shooting method which varies Ca* for 
constant initial conditions (that is, h = h' = h" = 0). We show in Figure the evolution 
of the critical capilary number with 8* , together with the diagram showing the different 
regimes described above. Notice that since 8* is a rescaled parameter, we have been able 
to investigate a large range of values, up to 8* « 10 7 . 

We did not find any reason for the disappearance of the contact line solution to coincide 
with the appearance of the LLD film solution as Ca* is varied. It may well be possible 
that, as W has already intersected the /i"-axis, the LLD film trajectory rolls up around 
the fixed point without h ever becoming negative. Some numerical simulations give us 
confidence that hysteresis indeed occurs (that is Ca* > Ca* 2 ): a slight hysteresis may 
indeed be observed in Figure HO Again, the present study is limited to LLD films of null 
flux, and other solutions may exist for the same parameters values. Thus, the hysteresis 
here observed can only describe a reduced part of the solutions set. 

When the capillary number is decreased from a supercritical value, the height of the 
stationary film hf decreases, and eventually the film thickness vanishes at some point 
Xmin (see Figure [5J ■ This point must be a minimum and in that case both h(x m i n ) and 
h'{x m in) vanish, so this film solution is also a contact line solution. This explains the 
change of the sign of h(x m i n ) observed at point A in Figure 

The asymptotic behavior of Ca* at large 8* has also been investigated (Figure 0). We 
find that the critical capillary number behaves asymptotically as a power law of the tilt 
angle, that fits to: 

Ca! ~ 0.3936 6* M " f 

e*->oo 

This suggest that for high 8* , Ca* behaves like 8* 3 / 2 



(4.3) 



5. Asymptotic results 

In the followin g section we describe a rough analytical approach, inspired from that 
of Eggcr which leads to the power law (|4.3(l for the capillary number obtained 

numerically in the previous section. 

5.1. Overview 

To determine the dependence of the critical capillary number with the angle, we need to 
better understand how the solution near the contact line connects with the free surface 
at infinity. We therefore seek to determine the matching between these two domains. 
This has actually been done when considering the classical plate withdrawal problem. It 
involves a matching between three regio ns: one near th e contact line, a capillary-viscous 
one and the gravity-capillary interface llEggersll2004l . However, two major differences 
arise in the present case compared to the usual problem. First, equation l|3.1|l is regular 
over the whole range h € [0;+oo[, whereas in the usual problem Navier slip condition 
l|2.3l) leads to a pressure divergence at the contact line. However, it has been shown (in the 
case of an advancing contact line) that the exact form of the slip law near the contact 
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Figure 6. Above: Second derivative of the film height against the capillary number, for 6* = 1. 

The solid line corresponds to h"(0) for the contact line solution. The dashed line represents the 

second derivative of h at the point a; m i n where in the case of a LLD film solution, the film is 

thinnest. Large dashes are used if h(x m i n ) > 0, and short dashes otherwise. 

Points A,B,C and D correspond respectively to the following values of Ca*: Ca* 2 , Ca* 1; 6* and 

0*/4. 

Below: Two LLD film solutions for 9* ~ 66.62 at different capillary numbers. The curve below 
correspond to point A, since h(xmi n ) becomes negative. 

line doe s not influence the matching proced ure with the far-field solution of the free 
surface (|Eggersll2004i iDussan fc Davislll974 . The second difference lies in the contact 
angle condition that we consider null instead of small but fin ite for solid plates. Such 
a condition is crucial as it can be seen from loggers! 12004^) where the solutions are 
expanded in powers of the small parameter Ca/#g (9 e being the static contact angle). 
Consequently, we cannot obtain a proper matching between the behavior in the contact 
line zone (cubic polynomial at leading order) with the famous logarithmic behavior in 
the capillary-viscous region: 

h ' [x) - [9Caln( ^k )]1/3 

However, we can bypass this difficulty by a slight change in the equations leading to 
a single approximation valid over the first two regions (contact line zone and capillary- 
viscous one). 

5.2. A two- zones matching 

The procedure hereafter presented is based on the assumption that the linear term in 
the right-hand term denominator of equation (|3.1|) is not of fundamental importance. In 
particular, the coefficient a was arbitrarily set to one at the beginning of this study, but 
numerical investigations have shown that it may be set to much different values, changing 
the results of only a few percent. For values of 1/a greater than 2/^/3, equation jHUJ 
becomes singular for negative values of h, which are physically meaningless. We therefore 
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Figure 7. Behaviour of W close to the largest fixed point . 0: no fixed point; +: positive fixed 
point; — : negative fixed point; a: spiraling trajectories; b: monotonic trajectories. Continuous 
lines delimit the various behaviors. The dashed line represents the critical capillary number 
above which the contact line solution disappear Ca*. Insert: asymptotic behavior of the critical 
capillary number for large 9* , (logarithmic scale). 

set from now on 1/q = 2/\/3 without any change in the equation properties, so that 
equation (j3.1J) reduces to 

h'"-h' + 9* = §^ . (5.1) 

(VV3 + 1) 2 V ' 

Close enough to the contact line, the film slope hi may be neglected due to the boundary 
condition l|3.1[) . In addition, close to the critical capillary number, Ca* ^> 9* (this is 
suggested by the assymptotic behavior H4.3|) 1. and equation (|5.1|) becomes 

h"'= (V 



(VV3 + 1) 2 ' 

solved analytically (as performed bv lDuffv & Wilson! lll997ln after the rescaling 

X = (3COV3 *' h{x) = ~ 1} - 

This rescaling leads to Tanner's problem: 

y'" = 4' y(°) = !> y'(°) = °- 

y 
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Its solution may be parametrized in terms of Airy functions Ai and Bi: 



r _ oi/3 Bi(s )Ai(s) - Bi(s)Ai(s ) = 1 

5 Bi'(« )Ai(s) - Bi(s)Ai'(s ) ' ^ 2 (Bi'(s )Ai(s) - Bi(s)Ai'(s )) 2 ' 

where sq is an integration constant, and s varies between consecutive solutions of 
equation 

Bi'(s )Ai(s) - Bi(s)Ai'(s ) = 0. (5.3) 

It was shown in Section f4.1.3l that. at the critical capillary number, h" vanishes (and so 
does y"). This property sets sq to zero and the range for s to [si, 0[, where Si « —1.98635 
is the largest solution to equation (|5.3(l . Matching solution (|5.2(l with the meniscus solu- 
tion should provide a condition on Ca* and 9*. At large £, the behavior of y ln is 

y m = a? + + c + O(-) 

where a,b and c have the following expressions (if we define z(s) = 7r(Bi(0)Ai(s) — 
Bi(.)Ai(0))) : 

Bi '<°> ^ 0.758947, ^afjS^U 1.12697, 



2 1 /3Bi( Sl )y ' ' VBi'(0)z'( Sl ) 

b 2 Bi'( Sl )s"( ai )-Bi"( Sl )s'( ai ) „ onfi?10 

C = ~a + W^W) 3 ' 

Matching y- m with the second-order Taylor expansion of for vanishing x leads to: 

9* - 0* p = 2a3 1/6 Ca* 2/3 , ^ = V3(c - 1) - 9* + 0* p , 6»* p = fe(3Ca*) 1/3 . 

When the critical capillary number tends to infinity, the first of the previous equations 
reads 

Cat - - * ,, A 9* z / 2 . (5.4) 

c e«-,oo (2a) 3 / 2 3 1 /4 v ; 

This fits remarkably well with the numerical estimation (|4.3() . since l/((2a) 3 / 2 3 1 / 4 ) w 
0.40631 (the above numerical fit giving 0.3936). The matching is compared to numerical 
results in Figure |H1 for a high value of 9* , and we observe a remarkable agreement with a 
reasonably large overlap region. 

Moreover, this matching procedure provides the following law for the apparent con- 
tact angle, at the critical speed: 9 ap — ^(SCac) 1 / 3 (notice that the rescaling term A 
disappears) . This is reminiscent of the famous Tanne r law which relates the contact line 
velocity to the apparent contact ang le l)Tannerlll979]) . 



6. Discussion and conclusions 

In this work, a continuum model of the forced dewetting on a porous material has been 
presented. In the framework of lubrication, a n ordinary non- linear differential equation 
was derived, close to the one investigated bv Irlockiria (|200 ll) . Even if the (microsopic) 
contact angle is assumed to vanish, a stationary contact line is found to exist for low 
dewetting velocity. Moreover, a transition between this steady contact line and the deposit 
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Figure 8. Comparison between numerical results (solid line) and the matching presented in this 
paper (long and short dashes), for 8* w 1886.19 and Ca* ps Ca* » 32871.1, in logarithmic scales. 
Long dahes: analytical solution to Tanner's problem; short dashes: capillary-gravity meniscus. 



of a LLD-film must occur, since there is a critical capillary number above which no contact 
line solution can exist. 

In the present study, the apparent contact angle (as denned by lEggersI ( 20(14 at) , that 



is #* p A in the present notations), does not vanish at the critica l capillary num ber (see 
Section [SJl . This behavior is different from the one obtained by lEggersI l|20044 . where 
the apparent contact angle was found to vanish at the transition from a contact line to a 
LLD film in the case of finite microscopic contac t angle d F . This discrepancy is surprising 
since one would expect the problem studied by lEggersI ( 2004ol) to correspond with the 



present one when 9 e — > and A — > 0. However, in the latter limit the rescaling by A 
used here can not hold. Also, as previously mentioned, the expansion in powers of Ca/ Q\ 
performed bv lEggersI l)2004a|) becomes ill-defined. 

Different c onclusions can be d rawn from our results regarding the erosion experiment 
performed bv lDaerr et all |2003). First, the existence diagram of the contact line can be 
drawn using the experimental values of the physical parameters. In the present theory, the 
permeability k of the porous material is crucial, as is the characteristic slip length at the 
solid- liquid interface. The value of this parameter may depend strongly on the compaction 
of the granul ar material ( s ay bet ween 10 -12 m 2 and 900 . 10~ 12 m 2 , respectively the value 
measured bv lDaerr et all l|2003t) and the square of the grain size). Figure El presents the 
critical velocities obtained for these two extremal values of the permeability. For the 
lowest permeability, and down to the smallest withdrawal velocities of the erodible plate, 
no contact line can exist. On the other hand, when choosing the largest permeability, the 
critical speed line is of order of those of the experiement. 

The flow acts on the granular medium mainly through the bottom shear r ate r = 
du/dy, which is known to trigger the erosion process (see lCharru et all l)2004|) ). From 
this shear rate, we can define the Shields number S, which compares the viscous force 
applied to the grains by the flow, to gravity force: 

a 7 n 



{p g -p)gd : 

where p g is the density of the grains (this expression stands only for small angles). 
Though the erosion prossess on a granular bed results from discrete and complex phe- 
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Figure 9. Critical speed for two values of the permability k. For each value of the permeabiblity, 
a contact line can exist only above the corresponding solid curve. The dashed line represents 
constant Shields nu mbers S c . The va lues of the physical parameters used in this model are 
those of Daer et al. (iDaerr et ffll.ll2003T) : g = 9.81 m.s -2 , p = 1000 kg.m~ 3 , p g = 2750 kg.m~ 3 , 
n — 10~ 3 kg.m _1 .s _1 and 7 S = 0.07 N.m -1 . The velocity and inclination ranges are those of the 
experiment. 



nomena, classical models assume that it starts at a threshold value of the shear rate, at 
which a critical Shields num ber S r is defined (Eh arm et all\2(}04. A typical val ue for 
S c is 0.05 (see among others iFredsoe fc Deigaardl l)l992l l h but IDaerr et all l)2003() used 
S c = 0.12 to fit their data. In addition to this large range of possible values, note that 
the critical Shields number is a function of the slope of the bottom: the more inclined 
it is, the easier it is for the flow to lift grains, thus the tilt reduces the value of S c . The 
shear rate may be deduced from our model, as a function of the dimensionless height of 
the film h: 

Uh 

Vk(h 2 /3 + h+l)' 

This expression admits a maximum value r max = U/(-\/k{\ + 2/V3)) for h = a/3, that 
is necessarily reached in the case of a contact line, since h stretches from zero to infinity. 
If the Shields number is assumed to be independent of the tilt angle of the plate (to first 
order), S = S c defines a vertical line in FigureEl (represented only for k = 900. 10 -12 m 2 ). 
On the left of such a line, no erosion should occur since the Shields numb er is smaller 
than t he critical value, thus to account for the erosion patterns observed bv lDaerr et all 
( 2003) at small velocities, a low value of S c is required. 

Now, if a LLD film covers the plate, h admit a minimum value h m i n which is found 
numerically to increase with the plate velocity. Once h m i n is larger than h = \/3, the 
maximum value of the shear rate becomes smaller than r max , and could decrease as the 
capillary number increases. Thus, the transition from contact line to LLD film could 
induce a strong change in the stress regime. If the physical parameters are choosen 
in the experimental range, no stress jump is observed numerically at the transition, 
but two assumptions should be relaxed in order to evaluate precisely the bottom stress 
close to the critical capillary number: the permanent regime, and the constraint of null 
outflow (Q = 0). In particular, if a negative outflow appears in reality (that is, if water 
is withdrawn from the bath), the shear rate should be reduced. Such a sharp stress 
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variation could provi de an explanation for the transition between the different erosion 
patterns observed bv lDaerr et all l|2003h . 

At the transition from contact line to LLD film, transient regimes should not be 
ignored. They have been studied in the litterature for non-vanishing contact angles 
( Hockina: l200lh . and future studies will aim at understanding the case of null contact 
angle, which has been shown here to be quite different. 

It is our pleasure to thank Daniel Lhuillier, Pierre- Yves Lagree, Eric Clement, Florent 
Malloggi and Jens Eggers for stimulating discussion. 



Appendix A. Pressure divergence at the contact line 

We aim here to demonstrate briefly that the classical Navier slip condition, leading to 
Equation J23J), is not sufficient to eliminate all the singularities at the contact-line, even 
for a non- vanishing microscopic contact angle (9 P > 0). Indeed, the first order expansion 
of Equation (E3J) is l|Eggerdl2004at lHockind ll983) 



n „ 3Ca ( , (xB e \\ 



where h is dimensionless but not rescaled (that is, h is the height of the water surface 
divied by the capillary length). In the lubrication approximation that we have used 
throughout, the pressure at the plate reads 



P\ y =o = pgcos(8)h--fh", 

and since h"(x) ~ — 3Ca/(0giz;) the pressure diverges at the contact line. In the present 
paper, due to the permeability of the porous plate, the expansion of h near the contact- 
line is a third order polynomial (h ~ h"(0)x 2 /2 + (Ca* — 9*)x 3 /6), and thus the pressure 
does not diverge. 



Appendix B. Derivation of the fundamental equation 

In the following, we aim to derive Equation (|2.(i|) from the two-dimensional Navier- 
Stokes equation, x and y refer to the axes of Figure ^ I n the frame of the lubrication 
approximation, and assuming both a permanent regime and small Reynolds number, 
momentum conservation reads 



1 Ov On 
— 7 f+gsm(9)+v— [ = (Bl) 
p ox oy z 

X dp 

— "o — gcos(6) = 0, 
pdy 

where u and v stand for the water velocity components respectively parallel and per- 
pendicular to the plate. The second equation may be integrated to give 



p = pgcos(8)(h -y)+p L 

where p^ is the pressure due to surface tension. Now, if 9 is small enough, the slope h! 
of the free surface should remain reasonably small, so that h" approximates its curvature, 
and pl ~ —"fh" (and, similarly, sin(0) w 9 and cos(9) « 1 at first order). The boundary 
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conditions on u at the bottom and the top of the film are (see Equations J23J, <|2.4|) and 

G3) 



Vk du 
U ^ + U= —d-y 



k I dp 

v =o v\dx 



y=o 



du 
dy 



= 0, 



y=h(x) 



thus integrating Equation (|B 1|) we obtain 



1 (\ dp 
v \p dx 



The mass flux is 



udy = Q. 



V. 



(B2) 



In steady state, mass conservation imposes that Q is a constant which vanishes in the 
case of a contact line. This condition provides the non-linear equation studied throughout 
this paper: Q = reads 



7 



ti" - gti 



uU 



P h 2 /3 + hy/k/a + k. 

which, after rescaling x and h by the capillary length l c , reduces to Equation (|2.6() . 
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